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ABSTRACT 

We study the Spitzer Infrared Array Camera (IRAC) mid-infrared (rest-frame optical) fluxes of 
14 newly WFC3/IR-detected z ~ 7 z 85 o— dropout galaxies and 5 z ~ 8 Yio5-dropout galaxies. The 
WFC3/IR depth and spatial resolution allow accurate removal of contaminating foreground light, 
enabling reliable flux measurements at 3.6/im and 4.5/xm. None of the galaxies are detected to [3.6] w 
26.9 (AB, 2a), but a stacking analysis reveals a robust detection for the zsso-dropouts and an upper 
limit for the Y105— dropouts. We construct average broadband spectral energy distributions using the 
stacked ACS, WFC3, and IRAC fluxes and fit stellar population synthesis models to derive mean 
rcdshifts, stellar masses, and ages. For the 0g5O — dropouts, we find z = 6.9j^J'i, (U — V) res t ~ 0.4, 
reddening Ay = 0, stellar mass < M* >= 1.2^0 6 x 10 9 -^0 (Salpeter initial mass function). The best- 
fit ages ~ 300 Myr, M/Ly « 0.2, and SSFR ~ 1.7 Gyr -1 are similar to values reported for luminous 
z rsj 7 galaxies, indicating the galaxies arc smaller but not much younger. The sub— L* galaxies 
observed here contribute significantly to the stellar mass density and under favorable conditions may 
have provided enough photons for sustained reionization at 7 < z < 11. In contrast, the z = 8.3^q 2 
Yios-dropouts have stellar masses that are uncertain by 1.5 dex due to the near-complete reliance 
on far-UV data. Adopting the 2er upper limit on the M/L(z = 8), the stellar mass density to 



Mijy AB < -18 declines from p*(z = 7) = 3.7 

-6 



+1.4 



X 10 6 Mr. 



m<3 Mpc~ A to p*(z = 8) < 8 x 10 5 M Q Mpc" d , 
following oc (1 + z)~° over 3 < z < 8. Lower masses at z = 8 would signify more dramatic evolution, 
which can be established with deeper IRAC observations, long before the arrival of the James Webb 
Space Telescope. 

Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

The reionization epoch represents the latest obser- 
vational frontier of galaxy formation. Little is known 
about the galaxy properties during this time, includ- 
ing feedback, metal production, and their contribu- 



tion to reionization. Extensive studies have been done 
of galaxies select ed by the dropout t echnique o ut to 
z = 6 fsee e.g.. IStanwav et al J 120031 lYan et all 120041: 
iBouwens et"alTl2006l iMcLure et all l2009af l. Pushing 
these studies to z > 7 has proven extraordinarily 
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rently known (|Bouwens et al.||2008 


: lOesch et al.ll2009at 


Castellano et al.ll2009l: lOuchi et all 


20091. R. J. Bouwens 



et al. in preparation). 

The arrival of WFC3/IR aboard Hubble Space Tele- 
scope (HST) heralds a dramatic improvement in our abil- 
ity to survey the reionization era, allowing us to iden- 
tify the predominant low-luminosity galaxies at z > 7 
by their rcdshiftcd UV light. Whereas HST remains 
crucial for identifying the galaxies, longer- wavelength 
data are necessary for constraining the stellar masses 
and ages. Access to the rest-frame optical is offered by 
the Infrared Array Camera fIRAC: iFazio et al.ll2004D on 
Spitzer, which by itself has permitted estima tes of stellar 
masses for large numbers of z ~ 5 — 6 sou rces (jEvles et al.l 
[20051 lYan et al.ll2006t iStark et al.ll2009h. and eve n a few 
at z ~ 7 (|Egami et al J 120051 : lLabbe et all 12001 ). One 
surprising early finding was the number of quite massive 
~ 1O 1O M0 galaxies with appreciable ages (200 — 300Myr) 
suggesting the galaxies formed substantial amounts of 
stars at even earlier time s , well into the e poch of reion- 
ization (jStark et al.H2007HYan e\~^MM \. 

In this Letter, we extend the stellar mass and age 
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Fig. 1. — Stacked images of z ~ 7 2g5o-dropout and 2 ~ 8 
ii05 -dropout sources in the HUDF. The stacked zg5o-dropouts 
in the [3.6]— band (top left) shows an unambiguous detection 
(S/N = 5.1). The median stack also shows a strong signal, indi- 
cating that it is not dominated by outliers. The S/N is calculated 
by bootstrap resampling. The zgso-dropouts stack is marginally 
detected (S/N = 2.6) in [4.5] (middle row). The Yios-dropouts 
stack (bottom row) shows no source in the center. Panels are in 
inverted grayscale and 9''9 X 9'.'9 on a side. 

estimates to lower luminosities (0.06 L*— 3 ) and higher 
rcdshift z w 8 by analyzing the IRAC mid-IR fluxes of 
the robust zsR n-dropout and Vms- dropo ut galaxy sam- 
ples fo und by lOesch et alJ (|2009bf ) and iBouwens et alJ 
(|2009aD in the newly acquired WFC3/IR data (GO- 
11563, PI Illingwort h) over the Hubble Ultra Deep 
Field ( HUDF; see also lBunker etail [20091 : iMcLure et all 
I2009H) . We adopt an fi M = 0.3, fl A = 0.7, cosmology 
with H Q = 70 km s^M pc" 1 . Magnitudes are in the 
AB photometric system (|Oke fc Gunnl fl983V Imaging 
depths refer to la AB point source total magnitude, un- 
less stated otherwise. 

2. OBSERVATIONS AND PHOTOMETRY 

The recent z = 7 2:350— dropout and z = 8 
^105 — dropout samples were selected from exceptionally 
deep WFC3/IR imaging, to « 29 AB magnitude 5a in 
Yio5,./i25,-Hi60, allowing the study of fainter sources than 
the N ICM OS based samples pr esented in lLabbe et alJ 
(|2006t) and lGonzalez et all (|2009f) . We discuss briefly the 
IRAC photometry o f 14 of the 16 zs5o~ dropouts found by 
lOesch et al.l (l2009b|) and the 5 Yins — dropouts re ported 
by IBouwens et al. ( 20093) (see lLabbe et al.ll2006l for the 



two Zg^o — dropouts found earlier with NICMOS). 

Spitzer data over this field were obtained from the 
Great Observatories Origins Deep Survey (GOODS; M. 
Dickinson et al. in preparation), which observed the 
HUDF with IRAC in two epochs of « 23.3 hr each 11 . The 
IRAC depths in the 3.6 and 4.5^m bands are 27.7 and 
27.2 (AB). We supplement the observations with deep 
if s -band data from the Very Large Telescope and Mag- 
ellan to K s ,ab = 27.4 (I. Labbe et al. in preparation). 

A critical aspect is obtaining reliable IRAC fluxes of 
the candidates. This is challenging because of extended 
wings of the IRAC point-spread function (PSF), caus- 
ing flux contamination fr om nearby fo r egrou nd sources. 
We use the technique of lLabbe et al.l (|2006f ) to fit and 
subtract the foreground sources using the flux profiles in 
the ultradeep WFC3 /IR detection images as priors (see 
Gonzalez et alj|2009l for a more complete description and 



Wuvts et al.ll2007l for illustrative examples). 



After cleaning the IRAC images, we first performed 
conventional aperture photometry in 3.6 and A.5fim 
bands in 2'.'5 diameter apertures. We define the pho- 
tometric error as the quadratic sum of random fluctua- 
tions in the aperture (determine d from empty aper tures 
on the nearby background; e.g.. lLabbe et al.ll2003D and 
the fit uncertainties to each individual neighbor. We per- 
formed the photometry independently on the first and 
second epoch IRAC data, which are rotated by 180 deg 
and therefore have different PSFs and contamination (see 
lLabbe et al.l 12001 ). Two of the 14 zsso-dropouts and 1 
Yios-dropout were too close to bright sources to be suc- 
cessfully measured and were removed from the sample 
(later we assume average M/L for these galaxies). Of 
the remaining galaxies none were detected (> 2a) in the 
epochs individually or combined (to [3.6] > 26.9, 2a to- 
tal). 

3. STACKING 

Owing to the extreme faintness of the zs50 — dropout 
and Fi 05 -dropout samples over the WFC3 HUDF09 
field, we derive average properties by registering and 
stacking the confusion corrected IRAC maps. The re- 
sulting stacked image for the 12 zsso-dropouts is shown 
in Figure [TJ featuring an unambiguous source coincident 
with the ^850— dropout location, with an average mag- 
nitude of [3.6] = 27.2 (AB, measured in a 2'.'5 diameter 
aperture and corrected to total) . To evaluate robustness 
of the detection we also create a median stack, which 
is resistant against outliers, to ensure that the signal is 
not coming from only a few sources. We derive uncer- 
tainties by bootstrap resampling the images in the stack 
200 times, finding signal-to-noise ratio S/N sa 5.1. The 
implied la rms in a 2'/5 aperture is [3.6] = 29.0 (total). 
The same procedure was carried out in the [4.5]— band 
showing a fainter source (SNR = 2.6). 

Detecting the Y105— dropouts is even more challenging. 
Not only are the Yios-dropouts in our samples fainter, 
there are also fewer to stack (4 versus 12), and they 
might be bluer due t o younger ages, lower metallicities, 
or other effects (e.g.. IBouwens et al.ll2009bh . The result- 
ing stacked image for the 4 is shown in the lower panels 

11 This paper uses data release DR3 of epoch 1 
and data release DR2 of epoch 2, ava ilable from 
http : //data. spitzer . caltech. edu/popular /goods/ 
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Fig. 2. — Average broadband ACS, WFC3/IR, and IRAC SEDs of the z ~ 7 and z ~ 8 galaxies. The left panel shows the average z ~ 7 
2850-dropout SED and the best-fit stellar population models. The blue far-UV slope and red Hieo — [3.6] color indicates a modest Balmer 
break, expected for evolved stellar populations (> lOOMyr). Excluding the IRAC measure ments leads to 1.5 dex s maller mass (dotted line). 
These faint galaxies have similar M/L and age as the luminous z ~ 7 sources presented by Gonzalez et al. (2009). The inset presents part 
of the SED around 1.4/im (~ 1700A rest-frame) show ing that the best-fit s ub-Solar metallicity models reproduce the quite blue far-UV 
slope better than do the best-fit Solar models (see, e.g.. IBouwens et~al . 2009b for more details). Differences between the best-fit parameters 
of various stellar population models at fixed metallicity are small. The Yios-dropouts (right panel) are undetected at [3.6]. The fit relies 
on the far-UV resulting in uncertain ages and masses. The gray dotted line shows the maximum M/L allowed by the fit (95% confidence). 
Clearly, moderately deeper IRAC would put stronger constraints on the allowed models at z = 8. Upper limits are 2a. Note, the models 
were fitted to the fluxes of the z ~ 7 and z ~ 8 galaxies, not the upper limits. 



of Figure [TJ showing no detection to a limit [3.6] > 27.5 
(2cr, total). The formal measurement is S/N = 0.9. The 
photometry is presented in Table 1. 

The average SED shape of the faint z ~ 7 
Z S50 —dropout galaxies, shown in Figure [2] is re- 
markably similar to that of the more luminous and 
massiv e NICMOS sample presented by iGonzalez et al.1 
20091). The rest-fram e far-UV slope is extremely blue 
Bouwens et "aT]|2009bl find (3 = -3.0 ± 0.2), suggesting 
low ages or sub-Solar metallicities. However, note the 
relatively red Hi eo — [3.6] = 0. 71^25 c °l° r > suggesting the 
presence of a modest Balmer break between the bands, 
indicative of evolved stellar populations. 

4. STELLAR POPULATIONS 

We use standard techniques to derive stellar masses 
and redshifts by fitting stellar populations synthesis mod- 
els to the average i, z, Y, J, H, K, 3.6, 4.5-band flux den- 
sities of the zrrq— dropou ts and 5*105 — dropouts. We use 
iBruzual fc Chariot! ( 
lar metallicity and a 
(IMF) between 0.1 



20031. BC03) m odels and assume So- 



Salpeterl (|1955| ) initial mass function 
and 100 Mq. This choice enables 
straightforward comparison with previous results, but 
we also review the effect of different assumptions and 
models. Furthermore, we assume constant star forma- 
tion (CSF) as opposed to exponentially declining models, 
motivated by the lack of evolution of the star formation 
rate (SFR) per unit mass (spec i fic SFR,(SSFR)) betwee n 
z = 3 and 7 (IGonzalez et all 12001 iStark et all [20091). 
Up to Ay = 0.5 mag iCalzetti et al.l (|2000i ) starburst 
reddening is allowed, where we note that z > 5 galax- 
ies have b lue UV-continuum slopes and low dust extinc- 
tion (e.g.. iLehnert h Bremer 2001 iStanwav et alJl2005t 
IBouwens et alll2?M I2009H) . We use the x 2 -fitting code 
FAST (jKriek et al.ll2009¥ which provides best-fit param- 



eters and uncertainties using Monte Carlo simulations. 
We fit the average SED fluxes to models smoothed to a 
resolution of 100 A rest frame, corresponding to the ap- 
proximate width of the dropout selection windows. The 
modeling results are presented in Table 2 and overplotted 
in Figure^] 

The best-fit redshifts of the zs50~ dropouts and 



I105— dropouts are z = 6.9 



+0.1 
0.1 



and 



8.3 



+0.1 



re- 
spectively, i n agreement with expectations from simu - 
lations (see lOesch et al.l I2009bb IBouwens etHI l2009al ). 
The extremely blue rest-frame far-UV slope of the z ~ 7 
galaxies is a challenge to fit, requiring Ay = 0, and sub- 
Solar metallicity rather than Solar models (Xred = ^.9 
versus x re d = 6-4). Consistently, the relatively red 
(U — V) res t ~ 0.4 color favors evolved models (age™ = 



350 



+30 
-170 



Myr ). We caution that emission lines can con- 
tribute to the [3.6] flux, but for the redshift distribution 
and iii6o ~ [3.6] color of this sample we estimate the effect 
is only ~ 0.05 — 0.1 mag. 

Interestingly, iMarastonl (|2005l ) and Chariot and 
Bruzual (in preparation; CB07) models produce nearly 
identical results as the BC03 models. Derived stellar 
masses are quite robust with at most 0.1 dex systematic 
variation between the different model choices. The only 
important effe ct on the masses is IMF: a more reasonable 
iKroupal (|2001[ ) IMF reduces the stellar masses and SFRs 
by 0.2 dex, but does not change other parameters or the 
quality of fit. More detailed exploration of the stellar 
ages and star formati on histories (SFHs) is presented by 
IGonzalez etahl (|2009l ). 

The stellar masses and ages of the Y105— dropout galax- 

12 Following ILabbe et al] (I200BT ) and IGonzalez et al.l 120091) , we 
report SFH weighted age™ = t/2 (for CSF), where t is the time 
elapsed since the start of star formation. 
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Fig. 3. — Evolution of the integrated stellar mass density. The 
red circle shows the z ~ 7 mass densi ty, derived by m ultiplying the 
integrated UV-luminosity density of[Bouwcns ct al. (2009aj) with 
the mean M/L of the sample studied here. The other z = 3 — 7 
results are based on luminous samples [M-uv AB < —20) from 
the lit erature (e.g.. [Sta rk ct al. 2009, blue circles: IGonzalez et al.l 
1120091 ) . green circle), corrected by +0.75 dex to Mijv,AB < —18, 
as appropriate for the UV LF of Bouwcns ct al. (2008) and a con- 
stant M/L. The dashed line shows a oc (1 + z)~ e evolution. The 
z = 8 upper limit shows the 95% confidence interval on the M/L, 
which is marginally in agreement with the extrapolated evolution. 
The floating error bar indicates the expected cosmic variance for 
the z ~ 7 and z ~ 8 samples. 

ies are not well constrained. The low S/N of the IRAC 
stack causes the fit to be determined by the J125 and 
Hiqq fluxes, producing a formal best-fit with minimal 
allowed masses ~ 4 x 10 7 M Q and ages ~ fO Myr and 
~1.5 dex uncertainties. However, if we refit the z ~ 7 
stack without IRAC data we also find very low masses 
and ages and large uncertainties. This underscores the 
crucial importance of deep IRAC data and suggests that 
model fits relying exclusively on far-UV data should be 
regarded with caution. Even though the fits may suggest 
a strong evolution of the stellar populations, we can not 
exclude that galaxies at z ~ 8 have similar masses, ages, 
and M/Ls as faint z ~ 7 galaxies. 

5. DISCUSSION 

The reionization epoch represents the latest fron- 
tier of galaxy formation theory, which is now being 
probed to extreme depths with the amazingly efficient 
WFC3/IR camera aboard HST. Using the conser- 
vative zssn— dropou t a nd Yi or— dropout sam ples of 
lOesch et ail (|2009bf ) and lBouwens et all (|2009af ) we find 
Spitzer /IRAC capable of an equally remarkable feat: 
the stacked 3.6/zm detection of even the faintest z ps 6.9 
galaxies found by WFC3 /IR, providing direct proof that 
IRAC can probe much deeper than widely accepted. In 
addition we derive upper limits for the newly discovered 
z pa 8.3 galaxies, some of which may well be at higher 
redshif t than the currently most distant z pa 8.2 ±0.1 
GRB (jSalvaterra et all [20091 ; iTanvir et al] [20091 ). We 
now elaborate on various implications. 



1 ) Stellar populations at z ~ 7 

Spitzer/IRAC places two important constraints 
on the stellar populations. First, the combina- 
tion of the blue far-UV slope and the moderately 
red Hiqq — [3.6] = 0.7+ ; 25 color in the ultrafaint 
(-H160 ~ 27.9) z ~ 7 galaxies indicates the presence of 
a modest Balmer break, expected for evolved stellar 
populations. Very young stellar ages are ruled out (CSF 
age^ > 80 Myr at 95%). Interestingly, low-mctallicity 
models, which produce bluer far-UV colors at fixed 
age, provide significantly better fits than do Solar, but 
for both it is challenging to reproduce the extremely 
blu e observed far-UV slo pe (/3 = —3.0 ± 0.2 as found 
by iBouwens et all l2009bf ) . Second, the average stellar 
masses < M >^ 1 x 10 9 M & indicate appreciable 
mass-to-light ratios M/Li 50 o « 0.1 and M/Ly ~ 0.2. 
The ratios are comparabl e those of the more m assive 
~ 7 x 10 9 M Q sample of IGonzalez et al.l (|2009t l. sug- 
gesting that the M/L ratios do not strongly depend on 
luminosity or mass. We find no appreciable difference 
between the best-fit parameters of BC03, MA05, and 
CB07 models, although IMF variations are a potentially 
significant source of s ystematic uncertainty: a steeper 
high-mass slope (e.g., IWeidner fc Kroupal I2005D would 
lower the inferred stellar ag es, a higher chara cteristic 
mass (e.g., "bottom light", Ivan Dokkuml 120081 ) would 
lower the total stellar masses and M/Ls, and any 
evolution with time would complicate comparisons 
between SFR and stellar mass across epochs. 

2) Star formation histories 

SFHs are notoriously hard to constrain from obser- 
vations of individual galaxies. However, statistical 
constraints can be obtained by comparing the SSFRs 
as a function of time and mass, as different shapes of 
the SFH give rise to distinctly different evolution (e.g., 
lLabbe et al.l 120071 ) . One remarkable recent result is 
the observation that the SSFR at fi xed ~ 5 x 10 9 M^ 
does not evolve fro m z = 3 — 7 (jStark et ail 120091 : 
Gonzale z et aill2009| ). Here we benefit from the extreme 
depths of the WFC3/IR data to probe the lowest 
luminosities, suggesting that the SSFR at z ~ 7 is also 
not strongly dependent on luminosity or stellar mass 
over the range 1 — 7 x 10 9 M Q . This may imply that 
the SFR correlates with stellar mass with a logarithmic 
slope close to 1. Both results (with redshift and mass) 
qualitatively agree with numerical simulations for early 
galactic star formation (z > 4), which robustly predict 
that SFR is proportional to M*, with si milar SFHs 
in haloes of different masses (|Davel I2008D . and that 
the SFRs of individ ual galaxies are constantly rising 
(jFinlator et al.ll200l . 

3) Evolution of the stellar mass density to z = 7 — 8 
The low-luminosity galaxies probed in this paper are 
expected to contribute significantly to the stellar mass 
density owing to their substantia l M/L . Following 
the approach of IGonzalez et ail (|2009f ) . we derive 
integrated stellar mass densities at z = 7 — 8 by mul- 
tiplying the UV-luminosity densities of IBouwens et ail 
(|2009al) , integrated to M uv . ab = —18, with the 
mean M/L derived for the galaxy samples in this 
Letter. For the ^850 —dropouts at z ~ 7, this results 
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in p* = 3.7t\\i x 106 MoMpc" 3 . The evolution of 
the mass density is shown in Figure [3l compared to 
previous luminous samples, which were corrected to 
the same luminosity limit assuming a constant M/L. 
The stellar masses at z = 8 are uncertain by ~ 1.5 dex 
due to the weak constraints from IRAC, so instead of 
plotting the best-fit M/L we show the 95% confidence 
upper limit. The integrated mass density to a limit of 
M UV ^ A B < -18 then declines to < 8 x 10 5 M Mpc -3 . 
This is marginally consistent with the extrapolation of 
the tx (1 + z)~ e evolution over the range 3 < z < 7. The 
upper limit may be evidence for a more rapid evolution 
toward z = 8, but deeper IRAC data and larger samples 
are needed to establish this firmly. The uncertainties 
quoted above do not include cosmic variance, which is 
estimate d to be ~ 30% for the z ~ 7 sample and ~ 40% 
at z ~ 8 (|Trenti fc Stiavellil[200l . 

4) Reinionization since z ~ 11 
Combined, the results have profound implications for 
the ability of the universe to be reionized by photons 
from star formation - and to remain so for an extended 
period. Crucially, the large stellar mass density in 
z ~ 7 low-luminosity galaxies points to significant 
amounts of star formation at earlier times. We can 
therefore ask whether this corresponds to enough 
ionizing photons at higher redshift to keep the universe 
reionized. Let us assume the stellar mass density 
p* = 3.7 x 10 6 M Q Mpc" 3 was assembled in the 340 Myr 
between 7 < z < 11 , in the reionization era (WMAP5, 
iKomatsu et al .1 l2009h . Correcting by 10% to account for 
mass loss in stellar evolution, the possible sustained SFR 
density would then be psfr ~ 0.012 M yr _1 Mpc~ 3 , 
significantly higher than that has been der ived for 
luminous galaxies alone ([Gonzalez et alJ 120091 ) . How- 
ever, the primary uncertainty is the fiducial critical 
SFR density neede d for reionization pff^ oc C/f esc 
(jMadau et al.l 119991 ). which depends sensitively on the 
values for the H II IGM clumping factor 1 < C < 30 
and the fraction of ionizing photons 0.05 < j e ^ r . < 1 that 
leak unhindered out of the galaxies (e.g., lOuchi et alJ 
I2009D . The popular choice (C = 30,/ esc = 0.1) 
leads to < p^f >= 0.7 M yr _1 Mpc -3 averaged 
over 7 < z < 11, which is 60 times higher than that 
can be explained by the stellar masses. However, at 
early times the clumping factor may be as low as 
C = 3-6 (e.g.. iTrac fe Cenl l2007t IBolton fc Haehnelti 
120071: iPawlik et al.ll2009fh The escape fraction may be 



significantly higher as well: recent numerical simulations 
suggest U.*r. > 0-2 fwlsefc Cedl2009t IPawlik et al.1 120091: 
lYajima et alTl2009[ ) and iBouwens et al.l (|2009bf f note 
that the extremely blue far-UV slope of the galaxies 
studied here could also result from large escape fractions 
fesc ^ 0.3. The choice (C = 3, f esc = 0.5) leads to 
< Peri? >~ °- 014 M yr- 1 Mpc- 3 (over 7 < z < 11), 
which means photons from the low-luminosity star 
forming galaxies observed here are capable of causing 
sustained reionization since z — 11. 

The first Spitzer/IRAC detection of low-luminosity 
z = 7 galaxies and their red H — [3.6] ~ 0.7 colors provide 
important insights into the earliest phases of galaxy evo- 
lution, showing that these early galaxies have relatively 
high ~ 300 Myr ages and M/Ly ~ 0.2. Consequently, 
low-luminosity galaxies contribute significantly to the 
stellar mass density p*(z = 7) = 3.7±ii x 10 6 M Mpc~ 3 , 
suggesting that the SFR corresponding to that mass 
may have provided a substantial fraction of the ionizing 
photons that kept the universe reionized at earlier 
times. Future IRAC studies of larger numbers of high 
redshift galaxies are required to bolster the results in 
this paper and to further constrain the evolution of the 
mass density toward the highest redshifts z > 8. The 
Spitzer post-cryogenic phase ( "warm mission" ) provides 
a rare window of opportunity to observe these galaxies 
in a timely manner to the ~ 0.5 mag deeper limits 
needed, long before the arrival of the James Webb Space 
Telescope. 
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TABLE 1 

Stacked photometry of z ~ 7 and z ~ 8 dropout galaxies in the HUDF 





-B435 


V606 


*775 


^850 


*105 


J125 


#160 


K s 


[3.6] 


[4.5] 


Z85n-dropout SED 


1.4(0.7) 


0.4(0.4) 


-0.1(0.6) 


5.6(1.1) 


27.2(1.8) 


29.6(2.0) 


24.5(1.6) 


22.9(13.8) 


48.5(9.5) 


43.2(16.6) 


Yios-dropout SED 


-1.1(0.9) 


0.2(0.7) 


0.3(0.8) 


-0.8(1.4) 


5.0(1.4) 


20.9(1.4) 


18.8(1.3) 


15.1(20.7) 


16.9(18.3) 


20.5(25.0) 



Note. — The optical-to- near-IR photometry is measured in 0'.'4 diameter apertures and corrected to total. Spitzer/IUAC photometry is 
performed on confusion corrected maps in 2'. 5 diameter apertures and corrected to total. Note, we limit the signal-to-noise ratio to S/N 15 
or 0.07 mag to account for zeropoint uncertainties. Units are nanojy. 



TABLE 2 

Best-fit stellar population parameters 
for constant star forming models and a salpeter (1955) imf 



model 


z phot 


Z log age™ Ay log M* 


log SFR 


log SSFR X l ed 






(Z Q ) (yr) (mag) (M ) 


(M yr" 1 ) 




2850-Dropouts Solar Metallicity 



BC03 




1.00 


8 co+0.00 


00+ ' 08 
u ' uu -0.00 


q 15 +0.07 
a - ±J -0.28 


q 7 +0.07 
-0.02 


-»■ '»-0.00 


6.39 


CB07 


6 - 85 -0'.07 


1.00 


8 co+0.00 


00+ 05 
u - uu -o.oo 


q 15+0.05 
a - ld -0.30 


n o 7 +0.05 
-0.02 


8 78+0.24 
'8-0.00 


6.39 


MA05 


f. o 7 +0.07 


1.00 


S - 4t) -0.36 


00+ ' 08 
u - uu -o.oo 


9 10+ 016 


43+0-08 
u -^°-0.02 


-8.67±° f 2 


5.48 



2g50-Dropouts sub-Solar Metallicity 



BC03 6.88±g;^ 0.20 8.58±g;g^ 0.00±g;£g 9.09±g;^ 0.3lt.g;J? -8.78± ; ou 2.95 
CB07 6.88±2'St °- 20 8 - 55 to.3p °- 00 ±o:oo 9 - 07 ±o.29 °- 31 to.oi - 8 - 76 ±!!:o2 2 - 94 
MA05 6.88tg;g 7 0.05 8.55±g;g| O.OOtqoo 9-07tg.ll °- 31 -o.oi -8-76to.il 2 - 72 



Y105 -Dropouts Solar Metallicity 



BC03 


8 91 +0.13 


1.00 


6 70+ 1 ' 77 
D -' u -o.oo 


00+ C) ' :iS 

u - uu -o.oo 


7 57+I.47 
'■°'-0.03 


0.58«;f 7 


-6 99+ 00 
o.»a_ 1 70 


2.31 


CB07 


8 n, +0.14 


1.00 


6 70+ 1 ' 77 
D -' u -o.oo 


00+ ' 40 
u - uu -o.oo 


7 c-7+1.51 
'•°'-0.03 


0.58«;f 7 


-6.99+??° 


2.31 


MA05 


8 9 o+0.15 
8 -^-0.23 


1.00 


6 70+ 1 80 
°-' u -o.oo 


00+ ' 45 
u - uu -o.oo 


7 =7+1.55 
'•°'-0.03 


0.59«; 4 3 


-6.98t?; 7 ° 


2.20 



Yi05-Dropouts sub-Solar Metallicity 



BC03 


8.26±°ii 


0.20 


6 70+ 1 - 80 

D -' u -o.oo 


o oo+°- 50 
u - uu -o.oo 


7 KQ+1-45 
'• OM -0.02 


60+ - 50 


-6.99+?' 1 " 


1.57 


CB07 


8 - 2 6t '.25 


0.20 


6 70+ 1 77 

D -' u -o.oo 


o oo+°- 50 
u - uu -o.oo 


'• oy -0.02 


60+ - 49 


-6.99l?:» 


1.57 


MA05 


8 - 2 6to-.25 


0.05 


6 70+ 1 77 
°- ' u -0.00 


o oo+ - 50 
u - uu -o.oo 


7 tro + 1.44 

'•°°-0.02 


60+ - 49 
U.DU_ 36 


_ 6 qq+0 00 
D - yM -1.70 


1.62 



NOTE. — Best-fit parameters and 68% confidence intervals as computed with FAST ( Krick ct al. 2009) 
for CSF models with a Salpeter (1955) IMF between 0.1 - 100 M Q . A Kroupa (2001) IMF would result 
in 0.20 dex lower stellar masses and SFRs, respectively, but other parameters would remain unchanged. 
The ^850 — dropout model fits favor low metallicitie s and high ages . The Y105 — dropout model fits are 
uncertain: deeper IRAC is needed. Note, following Gonzalez el al. (2009) we report the SFH averaged 
age w — t/2 (for CSF), where t is the time elapsed since the start of star formation 
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